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ABSTRACT 



Single and progressive -resolution coding algorithms for the 
compressian of 3-D polyhedral meshes are disclosed. In the 
single-resolution mode, the mesh topology (or connectivity) 
is needed by a constructive traversing api^ach applied to 
the dual graph of the original mesh while the mesh geometry 
is encoded by successive quantization and the bit-plane 
codii^ (achieved by context arithmetic coding). In the 
progressive-resolution mode, the mesh is represented by a 
coarse approximation (Le., the base mesh) and a sequence of 
refinements. Both the base mesh and the refinement opera- 
tions are entropy coded so that a series of mesh models of 
continuously varying resolutions can be constructed from 
the coded bit stream. Topological and geometrical data of a 
3-D mesh are encoded separately according to their impor- 
tance and theo integrated into a single bit stream. In 
decoding, the decoder finds firom the bit stream the most 
important information and gradually adds finer detailed 
information to provide a more complete 3-D graphic model. 
The decoder can stop at any point while giving a reasoiiable 
reconstruction of the original model. The disclosed algo- 
rithm was appUed to complicated 3-D meshes and adiieved 
a compression ratio of 20:1 while maintaining a good 
graphic quality. 
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3D MESH COMPRESSION AND CODING 
REFERENCE TO RELATED APPUCAnON 

This application is a conticuation-iD^art of provisional 
application Ser. No. 60A)73.087, filed Jan. 30, 1998, for 
Method and Apparatus for 3D Mesh Compresaon and 3D 
Mesh Coding 

BACKGROUND OF THE INVENTION 

1. Field of the Inventioo 

The present invention peitains generally to the field of 
computer graphics^ More particularly, the present invention 
pertains to coding of information representing an image of 
a three dimensional model 

2, Description of Prior Art 

Three dimensional ("3-Er) graphic models have become 
increa^ngly popular ance the advent of 3-D laser scanning 
systems and the boom of VRML (Wtual Reality Modeling 
Language) models. Laser scamnng systems, for example, 
routinely produce geometric models with hundreds of thou- 
sands of vertices and triangles, each of which may contain 
additional information such as color and normal. Highly 
detailed models are also commonly adopted in computer 
graphics. 

Three dimensional graphic models arc often represented 
as complex polyhedral meshes composed of two types of 
data, i.e., topobgical and geometrical data. Topological data 
provide coimectivity information among vertices (e.g., adja- 
cency of vertices, edges and faces) while geometrical 
attributes describe the position, normal, color, and applica- 
tion dependent information for each individual vertex. In 
terms of implementation, most 3-D graphic file formats 
consist of a list of polygois each of which is specified by its 
vertex indices and a vertex attribute list. The terms vertex 
and node are used interchangeably herein. 

Generally speaking, 3-D models are expensive to render, 
awkward to edit, and costly to transaiit through a network. 
Wder application of 3-D graphics potentially could be 
Umited due to these obstades. To reduce storage require- 
ments and transmission bandwidth it is desirable to com- 
press these models with lossy compression methods which 
keep the distortion within a tolerable level while maximizing 
data reduction. Another important consideration is to apply 
graphic coding in a progressive fashion to allow easier 
control of data such as progressive display, level-of-detail 
control and multi-scale editing. This functionality demands 
that the me^ be approximated with different resolutions, 
which is vital for real-time applications. To achieve this 
goal, it is required for the mesh to be reduced to a coarse 
approximatbn (i.e. the base mesfa) through a sequence of 
graphic simplifications. 

Simplification and compression of 3-D mesh data has 
been studied by quite a few researchers. Most early work 
focused on the simplification of graphic models. In W. J. 
Scfarocder, "Decimation of Triangle Me^ies,*" Computer 
Graphics Proceedings, Annual Conference Series, pp. 
65-70, ACM SIGGRAPH, July 1992, the author proposed a 
decimation algorithm that significantly reduced the number 
of polygons required to represent an object Turk, in "Re- 
tiling Polygon Surfaces," Computer Graphics Proceedings, 
Annual Conference Series, pp. 55-M, ACM SIGGRAPH, 
July 1992, presented an automatic method of creating sur- 
face models at several levels of detail &om an original 
polyhedral description of a given object. Hoppe et al., in 
"Mesh Optimization,** Computer Graphics Proceedings, 
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Annual Conference Series, pp. 19-26, ACM SIGGRAPH, 
August 1992, address the mesh optimization problem of 
approximating a given point set by using smaller number of 
vertices under certain topological constraints. 

^ Recent work has emphasized the compression of graphic 
models. Deering, in "Geometry Compression," Computer 
Graphics Proceedings, Annual Conference Series, pp. 
13-20, ACM SIGGRAPH, August 1995, discusses the con- 
cept of the generalized triangle mesh which compresses a 

10 triangle mesh structure. £ck et al. in "Muldresolution Analy- 
sis of Arbitrary Meshes," propose a wavelet transformatioD 
defined on an arbitrary domain to compress 3-D models of 
subdivi^n coimectivity. Taubin, in "Geometric Compres- 
sion Through Topological Surgery" Tedi. Rep. RC-20340, 

1^ IBM Watson Research Center, January 1996, presented a 
topological surgery algorithm which utilized two interleav- 
ing vertex aitd triangk trees to conipress a model. More 
recently, Cohen et aL, in ^Simplification Envelc^&,*' Com- 
puter Graphics Proceedic^ Armual Conference Series, pp. 

20 119-28, ACM SIGGRAPH, August 1996, introduced the 
concept of simplification envelopes so that a hierarchy of 
Icvcl-of detail approximations for a given polyhedral model 
could be generated automatically. Hoppe in "Progrcs«ve 
Meshes," Computer Graphics Proceedings, Annual Confer- 

25 ence Series, pp. 99-108, ACM SIGGRAPH, August 1996, 
proposed a progressive mesh compression algorithm that is 
applicable to arbitrary meshes. 

SUMMARY OF THE INVENTION 

30 

The present invention provides a new compression algo- 
rithm for 3-D meshes operating with both single and pro- 
gressive resolution modes. The single resolution mode com- 
presses the topological data through a constructive traversal 

33 and encodes geometrical data by local prediction. The 
progressive resolution mode represents the mesh by a base 
mesh and a sequence of refinement steps. Both the base 
mesh and the refinement sequence arc entropy coded into a 
single bit stream in such a way that, along the encoding 

4Q process, every output bit contributes to the reduction of 
coding distortion, and the contribution of bits decreases 
according to their order of position in the bit stream. At the 
receiver end the decoder can stop at any point while gen- 
erating a reconstruction of the original rnodel with the best 

45 rate distortion tradeoff. A series of models of continuously 
varying resolutions can thus be constructed firom the single 
bit stream. This property, often referred to as the embedding 
property, since the coding of a coarser model is embedded 
in the coding of a finer model, can be widely used in robust 
error control, progressive tran^nission and display, and 
level-of-detail control. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The other adjects, features, objects, and advantages of the 
present invention wiU be more fully understood and appre- 
ciated upon consideration of the following detailed descrip- 
tion of a preferred embodiment of the invention, presented 
in conjunction with the accompanying drawings, wherein: 

FIG. 1 is a graph illustrating the traversing order for each 
^ node for a given me^ in accordance with the invention. 

FIGS. 2{a) aiul 2(b) illustrates two link set configurations 
in accordance with the invention. 

FIG. 3 is a graph that illustrates the tracing of a merge link 
55 according to the invention. 

FIG. 4 is a graph ilhistrating the coding of a merge link 
according to the invention. 
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FIG. S illustrates a polyhedral me^ and its dual graph If every traversal step of a mesh is recorded, it can be 

according to one aspect of the invention. reconstructed by repeating the entire traversing procedure 

RGS. 6(fl), 6(b) and 6(c) are a graphical represeotatioD starting firom the initial node. Thus, during the traversal, a 

flhistrating recovery of a polygon according to the invention. link set is bookkept for each node of the mesh to record 

RGS. 7{a) and 7(A) are a graphical representation iUus- ^ relevant traversal results and assist the coding process. The 

trating vertex prediction according to the invention. configuration list cl(n) of node n stores aU links incident to 

FIGS. S(a) and m are a graphical representation illus- °«*f ° orders them in a counterclockwise order. The 

trating vert« spU? and edge coUapse according to the ^tol number of hnks m cl n) is the same ^ the valance of 

invecSon » ^ ^ ^o^^ ^ jjjg configuration list is cyclic m the sense that the 

^ . r ^ 1 **u 1 10 end of this list is implicitly connected to its beginning. The 

FIG. 9 IS a reprcsentahon of a 3-D model with one layer ^^^^^^ ^ ^nfiguralion list is the princq,le link. It 

of mdependent vertex ^ts. ^ . ^. ^ . serves as a starting point for tracing. Hie link list s(n> of 

FIG. 10 IS a graph ilhistrating the size distribution of 300 ^ ^ sequence of binary numbers which 

test models. denote whether a link is viatcd (labeled with 1) or not 

TIG. U is a graph ilhistrating the ratio of triangular (labeled with 0). Two Hnk configurations wdth respect to the 

meshes fiar 300 test models. central node arc shown in FIG. 2, where visited links arc 

FIG. 12 is a graph illustrating the coding performance depicted in solid lines and unvisitcd Knks in dashed lines, 

versus the number of vertices. We have cl(nHlo, li, I2, 1, 14, 15, Ig) and s(nHl. 0. 0. 1, 1, 

FIG. 13(fl) is a graph ilhistratii^ the performance of o, 1) for HG. 2(0), and cl(n)-(lo, Ij, 1^, 4) and s(n)-(l, 0, 

geometrical data coding with 10 bits/vertex. 20 1, O) for FIG. 2(b). The link set s(n) is modified whenever 

FIG. 13(ft) is a graph illustrating the performance of node n is involved in a traversal step, 

geometrical data coding with 20 bits/vertex. Coding Process 

FIG. 13(c) is a graph illustrating the performance of For leading node n^ visited in Step 2 of the constructive 

geometrical data coding with 30 bitsA^ertex. traversing algorithm, the modification of s(n,) is straight- 

FIG. 13(d) is a graph illustrating the performance of ^5 forward. Since new links arc visited in coimterclodcwise 

geometrical data coding with the rate-distortion tradeoff order, each of them is changed from 0 to 1 in the next tracing 

plot step sequentially. No extra coding bits arc required. 

FIG. 14(a) is a represenUtion of a model "dodge^ illus- each link Kn^ nj visited in Step 3 of the construc- 

trating compression at 10 bit^vertex. traversing algonthm, it can be a branch or a merger 

FIG. 14(b) is a representation of a model "dodge" illus- ^ depending on the type of n,. It is called a branch if n, has not 

trating compression at 15 bit^vertex. yf^^^'^ "'''^''^u^n P^^^^y- ^^Tf "^^^^1' 

„^ i*/v - . c ,u ^ 1 Abmary symbol bm(l) IS used to record the type of the hnk. 

FIG 14(c) IS a rcpresentation of the model dodge ^ ^J^^ ^ ^ yy ^ 

illustrating comprcssion at 20 blts^ertex. ^^^^ Furthe^ore, The valan^ v(nj is recorded and the 

HG. 14(d) IS a representation of the model dodge ^ initialized as (1, 0, . . . 1). For a merger e^n,, 

ilhistrating compression at 25 bits/vertex. n J, two additional pieces of information have to be collected 

DETAILED DESCRIPTION OF PREFERRED so that the traversal can be repeated to reconstruct the 

EMBODIMENTS original without ambiguity. They are the index id of node n^, 

Single-Resolution Mode and the position of e in s(n J 

The single resolution mode includes the coding of topo- 40 There are two ways to record the id of n^, i.e., a global 

logical and geometrical data. The topological data of a mesh indexing scheme, which is robust but costly, and a local (or 

are encoded losslessly. For geometric data, a local prediaioo relative) indexing scheme, which is efGcient but may fail 

and encoding of the residue is progressively performed to sometimes. To save coding bits, it is preferable to perform 

reach varying rate distortion performances. the local indexing scheme as much as possible. The global 

Coding of Topobgical Data--Constructive Traversal 45 indexing scheme is performed only when the bcal indexing 

A constructive traversing scheme is described below scheme fails. Let ^ be the set of links that have been 

where a mesh is constmcted by connecting a sequence of traversed already. The local indexing scheme attempts to 

traversal steps gradually. The traversing procedure begins locate n^ by tracing links contained by ^ in a clodcwise or 

with a single node and one link is traversed at a time. A counterclockwise fashion. The clockwise tracing procedure 

queue 0 of nodes is maintained during the traversal. It 50 is detailed bebw, where s(n)[i]issued denote the i** element 

records the order of eadi node to be visited. of s(n): 

The traversal is carried out based on 0 as follows: 1. Initialization: let n^ be n,, l^Knc* nj and set counter 

1. Initialize queue 0 with a node with index n. Pick one step to 0. Let 1^ be the ith element of cl(nc). 

of the links iocidem to n as its principle link. Both the node 2. Find the first j>i (in the cyclic sense) such that 

and the principle link can be selected arbitrarily. 55 s(n^)[j]=l and express cl(n^)[j] as (d^ nj. 

2. For the leading node n, in 0 traverse each of its 3. Increase step by 1, if i and j are not consecutive, 
unvisitcd links in a counterclockwise order starting torn its 4. Check whether n„«n^. If so, report step and terminate 
principle link. the tracing. Otherwise update node n^ to node n„. 

3. For each link l=<n^ nj in Step 2, if n. was not included This algorithm can be best explained widi an example as 
in 0 previously, add n, to © and set I as its principle link. 60 shown in FIG. 3, where visited links are depicted in solid 

4. After Steps 2 and 3 are completed for node n^ delete lines while unvisited links are depicted in da^ed lines. The 
node n, from 0. merger is shown in terms of the curved fine. Starting from 

5. If 0 is not empty, go to Step 2. Otherwise terminate the node n^, 11 the clockwise search takes five steps to locate the 
traversal. desired node n^ 12 based on links traversed before. The five 

A simple example is given in FIG. 1 to show die travers- 65 steps are illustrated with five arrow symbols. Note that each 

ing order of each node for a given mesh. This mesh step enck at either a node with a certain unfilled fink or the 

represents an open surface with two boles, (a) and (b). target node n^. The counterclockwise tracing is basically the 
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same as the clockwise tracing except that we change j>i to method to encode the prediction error. In an article by H. 

j<i (in the cyclic sense). For the same example, the backward Hoppe entitled "Progressive Meshes," Computer Graphics 

tracing lakes only one step to reach the target node. The Proceedings, Annual Conference Series, pp. 99-108, ACM 

(racing pattern which gives a smaller number of steps is SIGGRAPH, August 1998, the disclosure of which is hereby 

chosen for coding. There arc situations where the above 5 incorporated by reference, the author arranged vertices via a 

tracing algorithm may faD. For example, the local tracing sequence of vertex splits, and applied the delta prediction 

can be trapped in an infinite loop. Then the global iidcxing followed by the Hu^an coding. In Taubin et al., "Geo- 

scheme has to be applied. It is based on the observation that metric Compression Through Topological Surgery," Tech. 

the target node n^ has to be in queue 0 since (n^ n J is a Rep. RC-20340, IBM Watson Research Center, January 

merger. By viewing 6 as a cydic list, the ^obal method 10 1996, the authors defined a vertex tree structure so that the 

starts &om node n,, to search n^ in B top-down or bottom- position of a vertex is predicted by those of several ancestor 

up. The search is conducted one node at one step until it vertices along the tree. The prediction re^due is then trun- 

finds node n^. Again, the search direction which gives a cated and encoded with a Huffinan code. Choi et. al., in 

smaller number of steps is used and this step number is Results of Core Experiment M2/M3: Geometry Coding 

reported. To smnmarize, to ^jecify the id of the target node is Using PRVQ," Contribution Document M3148, MPEG4 

n^, we need a binary symbol for the search method (global San Jose Meeting, Fcbmary 1998, used the same vertex tree 

or local), a binary number for the search direction, and an and the prediction scheme but developed a multi-stage 

integer nuid>er for the number of search steps. vector quantization (VQ) method to encx)^ the prediction 

As far as the position of merger 1 in s(nj is concerned, it residue. The disclosure of the aforementioned articles is 

can be represented by an integer number counted from the 20 hereby incorporated by reference in its entirety, 

principle link in a forward or backward direction by con- Typically, there are three steps in a gieometry coding 

sidering visited links only. The actual counting direction is scheme: vertex ordering, data prediction and entropy coding, 

chosen to be consistent with the search direction for n^. In The first step arranges the list of vertices into a certain 

FIG. 4, an example is illustrated where node n, 13 has seven structure so that the local relationship among vertices can be 

links with the principle link Ig. Befiare the merge occurs, the 25 described and exploited more conveniently. The second step 

link list s(nj is (1, 0, 1, 1, 0, 0, 1, 1). There are only five utilizes such a structure to remove the redundancy of gco- 

visited links. Now, with link Ij as the merger, its portion is metrical data via prediction and produces a sequence of 

recorded as three when counted in a forward direction and prediction residues as a result The final step quantizes and 

two when counted in a backward direction. codes the residue according to a rate-distortion performance 

y^plication to Dual Graph 30 requirement. For simplicity, the focus herein is on the 

The dual graph G* 17 of a polyhedral mesh G 15 is compression of vertex positions. It is clear, however, that the 

constructed by converting a polygon in G into a node and same technique is applicable to other types of geometry data 

two nodes in die dual graph G* are linked if the correspond- such as surface normals and colors, 

ing polygons in the original graph G share the same edge. It The prediction of vertex positions is a causal process, ie., 

turns out that the dual graph G* is another polyhedral mesh. 35 a vertex has to be predicted by its preceding vertices. Any 

An example of a polyhedral mesh and its dual graph is prediction scheme should obey this rule so the decoder can 

shown in FIG. 5. There is a one-to-one correspondence decode the data properly. The definition of preceding verti- 

betwcen the original mesh G, 15 (dashed lines) and its dual ccs can, however, vary due to different proposed schemes. In 

graph G*, 17 (solid lines). Either one can be used to our proposed scheme, a vertex is predicted by some of its 

represent the topological structure of a 3-D graphic model 40 neighboring vertices. For a given vertex v, its neighboring 

uniquely. Since the above algorithm handles a general vertices are defined as those of polygons incident to v. 

polyhedral mesh it can be applied to the original me^ as However, not all of the neighboring vertices can be used for 

well as its dual graph. In the traversal of G the vertex valence prediction but only the preceding ones are eligible. The 

is distributed over a wider dynamic range which is typically preceding vertices of vertex v referred to herein are those 

from 3 to 8. However, most meshes primarily consist of 45 added to queue 0 before v. 

triangles and quadrangles so that the vertex valence is highly Vertex Prediction 

concentrated on 3 and 4. This property can be effectively There arc two approaches to predict the position of a 

ejqjloitcd by the entropy coder to make the proposed algo- vertex v. The first one regards v as one of the vertices of a 

rithm even more efficient. polygon and attempts to predict the position of v based on 

Coding of Geometry 50 those of preceding vertices of the same polygoa The second 

Geometry data specify attributes for vertices (or nodes) of one regards vertex v as the center of several vertices and 

a 3-D polyhedral mesh such as positions, normal, colors and predicts its position by the average of their positicms. The 

texture coordinates. They are functionally similar to pixel first approach is preferably used at first. If the first approach 

values of digital imagieA'ideo and sample values of digital is not applicable then the second approach is used. Since the 

audio. The major difference between digital audiovisual and 55 second approach is straightforward the following discussion 

3-D graphic media formats is that the former is usually provides more details on the first approach, 

defined on a regular spatial-temporal grid so that a transform By the first approach, a polygon is attempted to be rebuilt 

can be applied for energy compaction and redundancy based on a part of its vertices. It requires at least three 

removal while a 3-D graphic model is defined on an irregular vertices of a 3-D polygon to predict its remaining vertices, 

spatial grid Consequently its geometry data take the form of 60 For a plane determined by three points p^, p^ and p3 in the 

a list of sampled points on an arbitrary surface with no 3-D Euclidean ^ace, any point p in the plane can be 

particular order. Although geometry data bears a strong local represented by 
correlation it is difiEcult to find a simple transform to exploit 
such a correlation. 

The approach adopted currently for geometry compres- 65 where K^+X^+Xa^l. For a regular n-vertex polygon (p^pj - . 

sion is to use local praliction to remove data redundancy and . p^, if three vertices p,-, py and pjj arc known, then any other 

then apply a quantization scheme and a run-length coding vertex p/ can be recovered precisely by the above formula 
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where Kj, X3, depends only on i, j» k, I and n. To represent 
their dependent reUtionsbips of Xj, and X3, are denoted 
by a(n, i, j, k, 1), P(n, i, j,k, I) and ^(n, i, j k, 1), restively. 

The computation of a, and y can be done as follows. 
First, without toss of generality, it is supposed that the 
regular n-polygon is on tbc unit circle of the xy-plane such 
that 

PMf=oi (i^ sin (ifl)); 
Pf<<^ Oe). an Oa)X 
i>Xcos(l^suj(A^); 
Pr(a»(Za>,sm(yX 

where e=2^n. Then, we have 

{ctcos (cB)-pcos (jQ}vfcos (t:e><os (£6), 
{asin (»)=p5m (fe)n5in (ifi>-sin (©), 

The solution to the above system is: 



8in{(i' - 




-A)^ + sin((A 
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-*)^+siij((jt 




»in((i- 
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-m 


sin((/- 


+ sm((i - 




-m 



For an arbitrary n-vertex polygon (pip2 ... pj in the 3-D 
Euclidean space, the following approximation is adopted: 

In other words, p; is approximated by a linear combination 
of p,-, py and pj^ where coefficients arc set to those used to 
recover the regular n-vcrtex polygon of a ^>eci&c type. The 
recovered polygon can be mapped to any regular polygon by 
an affine transformation. FIG. 6 i^ows recoveries of a 
quadrangle, FIG. 6(a), pentagon, FIG. 6(6),and hexagon, 
FIG. 6(c), where the three preceding vertices 19, 21 and 23 
are depicted with black dots. These three vertices 19, 21 and 
23 are not required to be consecutive. If a polygon has more 
than three preceding vertices then any three of them 25, 27 
and 29 can be used to recover the polygon. To allow each 
preceding vertex to contribute equally, every possible com- 
bination of three preceding vertices is used to approximate 
the polygon. For a polygon with k preceding vertices, there 
are total of (k/3) different combinations. Each combination 
gives one prediction, and the polygon is finally recovered by 
the average of aU predictions. 

Since there are multiple polygons inddeat to a vertex, the 
position of this vertex can be predicted based on more than 
one individual polygon. In such a case, the final prediction 
is taken as the average of predictions produced by every 
applicable polygon. FIG. 7 shows the prediction of vertex v 
uiuicr two different circumstances, where its preceding ver- 
tices arc depicted with block spots. FIG. 7(a\ both the 
quadrangle 31 and the pentagon 35 can be recovered with 
predictions A, 37 and B, 39, respectively, and the average of 
A and B is taken as the prediction for v, 41. The above 
approach requires that each polygon contain at least three 
preceding vertices. This condition may not be met in some 
cases, such as FIG. 7(b). In such a case, the position of the 
vertex 43 is predicted by the average of all the preceding 
vertices. 
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For a model with only one connected component, the first 
vertex of queue Q has no ancestor. Its position can be coded 
either separately as overhead, or treated as the residue 
resulting from a null prediction. For a model with multiple 

5 connected components, the first vertex of tbc first compo- 
nent has no ancestor either. It is handled in the same way as 
in the single-component modeL However, the first vertex of 
any other component can be predicted by the position of the 
last vertex of the previous component Since the coder can 

10 freely dioose the coding order of components as well as the 
first vertex of each component, an optimai coding configu- 
ration to minimize prediction residues can be found. 
Successive Quantization and Bit-Plane Coding 

Most conventbnal graphic coding methods are devoted to 

IS the single resolution mode. Thai is, input cocffidents arc 
mapped onto a finite index set by a ^ngle quantization step. 
The index for each cocfiSdent is then converted to interme- 
diate symbols and encoded by an entropy coder. An embed- 
ded coding scheme as described herein may be used instead 

20 of the single resolution mode by adopting a different quan- 
tization procedure called succesave quantization. Succes- 
sive quantization uses a sequence of gradually refined quan- 
tization steps instead of a single quantization step. The 
quantization step size is refined so that coefficients are 

25 approximated by an increasir^ precision. For each quanti- 
zation step, a binary output is produced for each cocffident 
to record the quantization result All symbols generated via 
the same quantization step are then grouped together and 
entropy coded before the quantization of the next step. This 

30 is often known as bit-plane coding. Bit-plane coding can be 
effectively achieved by using the context arithmetic coder. 
For details of successive quantization, bit-plane coding and 
context arithmetic coding the reader is referred to Li et al., 
"Progressive Coding of 3D Graphic Models," Proceedings 

35 of IEEE, Vbl. 86, pp. 1052-63, June 1998, the disclosure of 
which is hereby incorporated by reference in its entirety. It 
is worthwhile to point out that even though the method 
presented above is called the single resolution method it still 
generates an embedded (or progressive) bit stream due to the 

40 embedded coding scheme applied to geometrical data. The 
term of "a single resolution" is adopted here to indicate that 
the topotogical structure (or connectivity) of a 3-D me^ 
remains the same throughout the coding process. 
Coding of Overhead 

45 To efficiently compress small as well as large-size models, 
the overhead bits preferable should be reduced to a mini- 
mum. First of all the initial quantization threshold Tq must 
be kept as the overhead. It is usually stored as a 32-bit 
floating point number, and can be reduced to as low as a 

50 16-bit floating point number to meet the low bit rate require- 
ment for the coding of models of small size. The position of 
the first vertex of the first component may be kept as the 
overhead also. It requires 3*32=96 bits. In an example of bw 
bit rate coding (10 bit/vertex) of a model of 20 vertices the 

55 total number of overhead bits is between 112 and 128 which 
is in fact more than one half of tbc entire bit budget (e.g., 200 
bits). 

To improve the situation, the position of the first vertex 
may be coded as a residue. In other wwds, it is successively 

GO quantized and encoded (in contrast with the codii^ of its 
whole value at once). There are two consequences of this 
chokx. First, the position value is usually much larger than 
prediaion residues. By treating the position of the first 
vertex as a residue, the magnitude of the initial threshold Tq 

65 has to be increased ^gnificantly and the first several quan- 
tization steps are solely used to encode the residue of the first 
vertex. For small size models, this choice is still more 
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efficient than storing the position as a floating-point number 
directly in the overhead For large size models this choice is 
less efficient and costs more bits. However, such a cost is 
negligible compared to the relatively large bit budget for 
large size models. Second, the first vertex serves as an 5 
anchor point for the entire mesh. The mesh could jump 
around as the position of the first vertex is refined gra^ally 
in the successive decoding stage. However, it does not 
introduce any visual distortion since the shape of the mesb 
is not deformed. lO 
Progressive-Resolution Mode 

Although the single-resolution compression method 
described above achieves a significant coding gain for 
topological and geometrical data reduction it may not be 
adequate in some cases. A simple analysis of the comprcs- 15 
sion performance is given below. For a closed triangular 
mc^ with o vertices and approximately 2n triangles, each 
vertex requires 32x3=96 bits to store the x, y and z coordi- 
nates as a Soating-point number aixi each triangle requires 
32x3"^ bits to store the three vertex indices as an integer 20 
number. If the mesh is represented in the binary format the 
number of bits required may be represented by the following 
relation: 

Rx964^2nx96=i2S8nfaits. 25 

The total number of bits can be 3~5 times larger, if the mesh 
is represented in ASCII format. When compressed, topo- 
logical (or connectivity) data requires approximately 2 bits 
per triangle while the position of a vertex takes about 20 bits ^ 
to generate a visually indistinguishable approximation. The 
total bit consumption is: 

/o(20+2nx 2*-24«bits 

Thus, the compression ratio is around 12:1. For a large 35 
model with tens of thousands of vertices, such a compres- 
sion ratio is still not high enough to meet some storage and 
transmission requirements. Furthermore, the single- 
resolution method attempts to reduce the number of bits for 
model representation, but makes no effort to simplify the 40 
graphic model so that its retxiering and editing speed can be 
enhanced. Real-time interactivity cannot be easily achieved 
for complicated single-resolution models even with today's 
high performance computers. For many applications, real- 
time rendering is frequently more important than the com- 45 
prcssion issue. To overcome these two problems the 
progressive-resolution representation of a 3-D mesh pro- 
vides a good solution. 

As described herein, the progressive-resolution mode is a 
3-D mesh that is gradually tran^ormed into a sequence of 50 
simpler me^es through graphic Amplification steps. The 
final and simplest mesh is often referred to as the base me^. 
To allow the conversion from one resolution to another, 
every simplification step has to be inevitable. This inverting 
process is called refinement. Thus^ the original mesh can be S5 
reconstructed by the base mesh and a sequence of refine- 
ments arranged in the reverse order of the simplifications. 
Any intermediate resuh in the process can be used as an 
approximation of the original mesh. The more refinement 
steps utilized, the closer the approximation is to the original 60 
mesh. This multi-resolution property is related to topologi- 
cal data coding. However, the multi-resolution technique 
can also be applied to geometrical data coding as presented 
in the section on Successive Quantization and Bit-Plane 
Coding. These two progressively coded bit streams may be 65 
integrated into a single embedded bit stream for 3-D mesb 
representation as discussed below. 
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Graphic Simplification 

The polygon is the primitive of an object represented by 
a polyhedral mesb. The cost of storing or transmitting a 
polyhedral mesh of n primitives is 0(n). The cost of ren- 
dering such a model is also 0(o) as discussed in Bar-Yehuda 
et al., "Tmae/Spaoc Tradeoffs for Polygon Mesh Rendering," 
ACM Transactions on Graphics, Vol. 15, pp. 141-52, April 
1996, and in Heckbert et aL, "Multiresohition Modeling for 
Fast Rendering,** Proceedings of Graphics Inter&ce '94, pp. 
43-50, Canadian Information Processing Society, May 
1994. The disclosure of each of these two articles is hereby 
incorporated by reference in its entirety. 

The primary goal of graphic simplification is to reduce the 
number cf primitives required to faithfully represent an 
object Several different algorithms have been proposed to 
serve this purpose. They can be roughly classified into three 
categories i.e., surface retiling, vertex decimation, and ver- 
tex clustering, as described below. 

With surface retiling, polygonal surfaces are triangulated 
with a new set of vertices to replace the original set The new 
set usually consists of fewer vertices than the original one 
while preserving its topology, in Turk, "Re-tiling Polygon 
Surfaces," Computer Graphics Proceedings, Annual Confer- 
ence Series, pp. 55-64, ACM SIGGRAPH, July 1992, it is 
suggested to randomly add points on the polygonal surface 
and then re-distribute those points by exerting a repelling 
force while removing old vertices graduaUy. Hoppe et al., in 
an article entitled, "Mesh Optimization," appearing in Com- 
puter Graphics Proceedings, Annual Conference Series, pp. 
19-26, ACM SIGGRAPH, August 1993, defined an energy 
function which is minimized to determine positions of new 
vertices. In Hinker et al., "Geometric Optimization,** Proc. 
Msualization, pp. 189-95, October 1993, coplanar and near 
coplanar polygons are merged into larger complex polygons 
and then re-triangulated into fewer simple polygons. Kalvin 
and Taylor in "Superfaces: Polygonal Mesh Simphfication 
with Bounded Error," IEEE Computer Graphics and 
implication, \bL 16, pp. 64-77, May 1996, developed a 
similar algorithm which allows additional control of the 
approximation error. In Cohen et al., "Simplification 
Envelopes," Computer Graphics Proceedings, Annual Con- 
ference Series, pp. 119-28, ACM SIGGRAPH, August 
1996, and in Varshney et al., "Generating Levels of Detail 
for Large-Scale Polygonal Models," T^h. Rep., Department 
of Computer Science, Duke University, 1995, the original 
polygonal surface is surrounded with two envelopes and 
then generated a Amplified surface within the volume 
enclosed by these two envelopes. The disclosure of each of 
the aforementioned articles is hereby incorporated by refer- 
ence. 

By using vertex decimation, Schroeder, in "Decimation of 
Triangle Meshes,** Coniputer Graphics Proceedings* Annual 
Conference Series, pp. 65-70, ACM SIGGRAPH, July 
1992, proposed to make multiple passes over an existing 
polygonal mesh and use local geometry and topology infor- 
mation to remove vertices which meet a distance or an angle 
criterion. The hole left by the vertex removal is patched by 
a local triangulating process. Soucy et al., in "Multiresohi- 
tion Surface Modeling Based on Hierarchical 
Triangulation," Computer Vision and Image Understanding, 
Vot 63., pp. 1-14, January 1996, describe a more sophisti- 
cated algorithm following the same framework. The disclo- 
sure of each of the aforementioned articles is hereby incor- 
porated by reference. 

The idea behind vertex clustering is that a detailed part of 
a model is represented by a set of spatially close vertices. A 
clustering algorithm is used to generate a hierarchy of 
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clusters and approzimatxons with different resohiticxns. For a 
discussion of this concept see Schaufler et al., "Generating 
Multiple Levels of Detail £roai Polygonal Geometry 
Models,'' Mrtual Environments, Eurographics Workshop, 
pp. 33-41, January 1995. Rossignac et aL, in "Modeling in 
Computer Graphics: Methods and Applications,'' pp. 
455-65, Springer- Vferlag, 1993, divided the bounding box of 
the original model into a grid. \^^thin each cell, vertices are 
chistered together into a single vertex, and the model surface 
is updated accordingly. He ct al., in "\foxel Based Object 
Simplification," Plroc. \^alization, pp. 296-303, 1995, 
adopted a signal-processing approach to sample the iiqmt 
object and used low-pass filtering to remove high fiequen- 
cies of the object \^^x split/edge colla^jsc is the most 
commonly used clustering technique. It removes an edge 
and contracts two end points into a new vertex. Tb determine 
the position of the new vertex, Gueziec, in an article entitled, 
"Surface Simplification Inside a Tolerance \felume,'' Tfech. 
Rep. RC-20440, IBM Watson Research Center, 1996, 
required the volume to be preserved after contraction. Hoppe 
in "Progressive Meshes," Computer Graphics Proceedings, 
Annual Conference Series, pp. 99-108, ACM SIGGRAPH, 
August 1996, and in "Progressive Simplicial Complexes," 
Computer Graphics Proceedings, Annual Conference Series, 
pp. 217-25, ACM SIGGRAPH, August 1997, used edge 
collapse to form the progressive representation of triangu- 
lated geometric models. Garland et a]., in "Surface Simpli- 
fication Using Quadric Error Metrics," Computer Graphics 
Proceeding?, Annual Conference Series, pp. 209-17, ACM 
SIGGRAPH, August 1997, developed a method called pair 
contraction which is capable of contracting an arbitrary pair 
of vertices. The disclosure of each of the aforementioned 
articles is hereby incorporated by reference. 

Although other techniques may be used, the vertex split/ 
edge collapse technique is preferably used for graphic 
simplification due to the following reasons. First, vertex ^ht 
naturally leads to a progressive representation of a model 
without any constraint. Second, an instance of vertex spht 
requires the least number of bits to encode among existing 
techniques. It can be coded very efiScicntly since the vertex 
split operation only involves two incident edges. Third, 
vertex sp]iX provides the finest granularity for simplification 
or refinement. Each operation only renders a small change of 
the model, and such a change can be realized continuously 
by a morphing process. This enables a smooth transition 
between different resolutions of a model. 
Vertex Split and Edge Collapse 

An edge collapse operation ecol({v„ v^, v^ v^}) removes 
ed^ (v„, v^ from the mesh and merges the two end points 
v^ 45 and 47 into a new vertex v^, 49 as shown in FIG. 
8. Edges previously incident to either v„ or Vj are then 
connected to vertex v^, IWo adjacent triangles (v^v^vj and 
(v^v^v^, FIG. 8(fl) vanish in this process. Each edge col- 
lapse decreases the number of vertices by one and the 
number of triangles by two. As a reverse process, a vertex 
split operation vsplit ({v;, v^., v^}) divides a vertex v^, 49 into 
two vertices v„, 45 and v^ 47. Consequently, edge (v„, v^ 
and two adjacent triangles (v,, VjV„) and (v^v^v^) are created. 
Each vertex split increases the number of vertices by one and 
the number of triangles by two. 

Hie edge collapse operation can be repeatedly applied to 
an original mesh denoted by M^ where one edge is 
removed at a time. The process continues until a simple base 
mesh M^^ is obtained. This process can be expressed as: 
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W— =W* -M^-l... -Ml ^Afo=W4«r 

" -* — ► -» 

^ The sequence of edge collapse operations should be chosen 
carefully so that each M,> approximates the original one with 
the best achievable quality. This problem will be addressed 
below. Each edge collapse operation can be inverted by a 
corresponding vertex split operation. Therefore the original 

10 mesh M^^ can be reconstructed based on the base mesh 
Mj^, and a sequence of vertex ^lit operations arranged in 
die reverse order of edge collapse operatic»is, i.e.: 

This forms a natural progressive representation of the 
original mesh. The base mesh Mq provides the coarsest 
approximation while mesh M„ is identical to the original 

20 one. Meshes M^, 0<i<n, incrementally converge to the 
original one as more vertex split operations are applied. 

Two edge collapse operations ecolCfv^^, v^^ Vj^ v^j}) 
and ecol({v2/, Vj^, Vj^}) are independent of each other if 
they do not ^are any common vertex. In other words, the 

25 intersection of sets {v^^ v^^, v^^ v^} and {v^/y v^^, v^J is 
empty. Two vertex ^lit operations are independent if their 
corre^ndit^ edge collapse operations arc independent. 
Based on this concept, a mesh simplification may be con- 
ducted through several consecutive layers where indepen- 

30 dent edge collapses are performed in the same layer. The 
independence constraint forces edge collapses ^read evenly 
over the surface of the model. Therefore, the mo^l is 
simplified across its entire surface in a balanced way. It 
ensures that the simplified model is of reasonably good 

35 quality. FIG. 9 illustrates one layer of independent edge 
collapses of model Triceratops, where each vertex split is 
plotted as a patch of two gray triangles. These patches are 
disjointed with each other. FIG. 9 also shows that two 
independent vertex splits are commutable. That is, they yield 

40 the same final result regardless of the order in which they are 
applied. For example, every patch in FIG. 9 can be collapsed 
freely without interfering with the collapses of other patches. 
A particular order of edge collapse can be chosen to give the 
best coding performance for intermediate meshes, as 

45 detailed below. Each edge collapse results in loss of infor- 
mation and introduces a certain amount of distortion. To 
optimize the rate-distortion performance, the edge with the 
least importance is collapsed at each simplification step. 
This process may be performed as described in Garland et 

50 al., "Surface Simplification Using Quadric Error Metrics," 
Computer Graphics Proceedings, Aimual Conference Series, 
pp. 209-17, ACM SIGGRAPH, August 1997, and Hoppe, 
"Progressive Meshes,"* Computer Graphics Proceedir^, 
Annual Conference Series, pp. 99-108, ACM SIGGRAPH, 

55 August 1996. The simplification step continues until either 
a certain quality criterion is met or any further simplification 
violates the topology of the model. To look from the vertex 
spilt viewpoint, one sbotild use the vertex split operation to 
recover the most important vertices first (such as the peak of 

60 a spike) and then less important vertices (such as vertices 
coplanar with their neighboring vertices). 
Coding Issues 

Base me^ and the sequence of vertex split opera- 
ticms {v^lit;, . . . , vsplit„} are coded in order to form a 
65 single embedded bit stream. The decoder decodes the bit 
stream and recovers M^^, vspKt^ . . . ,vsplit„ sequentially. 
From this single bit stream, the decoder can reconstmct the 
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model with continuous resolutioos, ranging from the coars- 
est approximation to the finest replica. Depending on 
application, the decoder can stop at any position of the bit 
stream to decode a mesh which approximates the original 
one to a certain degree. For example, a user may want to 
render a far-away model with a very coarse approximation. 
In this case, the decoder can stop right after it rebuilds the 
base mesh. On the other hand, if a user needs a close view 
of the model, the decoder may decode the entire bit stream. 
Furthermore, if a user is animating a zooming process by a 
sequence of frames, the decoder can incrementally decode 
the bit stream in accordance with the decreasing distance 
between the camera and the model to get the optimal 
rendering quality with the ininimal computational complex- 
ity. 

In comparison with onginal mesh base mesh 
is very simple. For example, a high quality mesh of the 
sphere may have tens of thousands of vertices and triangles. 
However, its base mesh is the simplest 3-D mesh, i.e. a 
tetrahedron. The base mesh can be coded as a regular mesh 
with the single-resolution mode. Only a few bits are required 
to encode the base mesh. 

The vertex split operations are applied in a layered 
fashion. The first layer starts from the base mesh, and every 
following-layer starts fi^om the mesh resulting from the 
previous layer. At the beginning of each layer i, all vertices 
are ordered as a ID list L.- and every vertex in the list can be 
split into two new vertices. List L,- is modified when a vertex 
split occurs. One binary flag splitFlag is allocated for each 
vertex. If vertex v indeed ^lits, splitFlag(v) takes the value 
of TRUE. Otherwise, it takes the value of FALSE. The 
coding of layer i starts from the first vertex of list L^-. One 
vertex is processed at a time. Let v^ denote the vertex under 
current procesang. If v^ does not split, splitFlag (v^) is set 
to FALSE and coded, the en£X)der gpes to the next vertex in 
Li, If v^ is indeed split by vsplit({v,, v^, v^), splitFlag (vj 
is set to TRUE and coded. Furthermore, the vertex id of V; 
and v^can be coded jointly. Suppose v^ is of valence n before 
splitting, and V; and v^ are two of its n neighboring vertices. 
There are n(n-l)/2 different patterns to choose two out of the 
n neighbors. It is possible to code v, and v^ jointly from a 
careful analysis of these patterns. 

After v^ splits into two vertices v„ and v^ v^ replaces v^ 
m hi while v^ is added to the end of L^. TTien the encoder 
visits the next vertex in L^. This process is repeated until all 
vertices initially listed in L, arc axled. New vertices created 
during the coding of layer i and appended to L,- are not 
visited at the same layer. At the end of coding of the layer 
i, L,- becomes L^^^ and then the coding of layer i+1 starts. 

The position of new vertex v^ for an edge collapse can be 
determined and recorded in two ways. The simplest way is 
to place v^ in the middle of the collapsed edge (v^ v^. A 
more sophisticated way is to place v^ in such a position that 
the distortion due to the edge collapse is minimized. For this 
case, the new position is generally not the middle point of 
edge (Vy, vj. The second approach provides a graphic model 
of better quality at the expense of a higher computational 
complexity. Furthermore, since the distortion incurred at any 
stage will inevitably propagate to the latter stage the quality 
of the mesh deteriorates at a greater rate by using the first 
approach. However, from the coding viewpoint, the first 
approach requires fewer bits. Two di^lacements v„-v^ and 
v^-v^ must be encoded to perform the vertex ^lil operation 
with the second approach while only one displacement is 
needed with the first approach due to the fact that v„-v^a 
v^-v^. Since the coding of displacements takes the lion share 
of bit consumption, the first approach reduces the bit 
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requirement significantly. Under the same bit budget, it can 
encode many more vertex ^ht operations than the second 
approadi. The tradeoff is that under the same bit budget the 
first approach etKXKles more vertex splits while each vertex 

5 Split operation contributes less to improve the quality. It is 
observed that the fir^ approach gives a bctto* overall per- 
formance. Additional, the first approach is computationally 
simpler. Under the first approach, for each vertex spht v^tit 
(V;, v^, V J, the displacement of position v^-v^-v^-v^ must 

10 be coded in order for the decoder to rebuild the model. First, 
all vertex displacements are collected into a sequence in the 
order of occurrence of vertex splits. Then, the se q uence is 
coded by successive quantization and embedded coding 
techniques as described in li et al., '"Progressive Coding of 

15 3D Graphic Models," Proceedings of IEEE, VbL 86, pp. 
1052-63. June 1998, the disclosure of which is hereby 
incorporated by reference in its entirety. 
Integration 

Two coding procedures were presented before, i.e., the 

20 coding of cotmectivity and geometry data. Each coding 
scheme produces its own bit stream. These two bit streams 
have to be multiplexed into a single bit stream. The quality 
of the reconstructed model depends on both the number of 
vertices and the precision of vertex position. The more a bit 

25 stream is decoded, the more similar the decoded model is to 
the original model. Decoding of either bit stream contributes 
to the reduction of compression distortion. The contribution 
of bits decreases according to their order in each bit stream. 
It is desirable that this property is preserved in the final bit 

30 Stream as well. This problem is similar to the merger of two 
arrays, each of which has been sorted according to a certain 
measurement individually into one array. In this case, the 
measurement is *"distortion reduction per bit.^ A rate- 
distortion model is built for each bit stream to study the 

35 average distortion reduction per bit so that these two bit 
streams can be nmltiplexed in a proper order. For more 
details on this process, the reader is directed to the afore- 
mentioned article by Li et ah, entitled, "Progressive Coding 
of 3D Graphic ModeU." 

40 Experimental Results 

A test data set consisting of 300 models was used to 
validate the coding approach disclosed herein. The size of 
these 300 meshes is well distributed. FIG. 10 illustrates the 
size of each mesh in terms of the number of vertices. Note 

45 that the x-axis is plotted in the log scale. The smallest size 
is four, which corresponds to a tetrahedron while the largest 
size is 54898 which corre^nds to a complicated gear set. 

As mentioned before, most meshes consist primarily of 
triangles and rectangles. The ratio of the number of triangles 

50 to the total number of polygons for each mesh is plotted in 
FIG. 11. As can be seen from FIG. 11, the triangle is the most 
common form of polygon, the rectangle is a distant second 
and polygons of higher degree are relatively few. This can be 
explained by the fact that a plane in the 3-D EucUdean space 

55 is determined by three points which form a triangle. To 
represent a 3-D graphic model by a polyhedral mesh, the 
surface of the model is actually approximated by pieccwisc 
plane segments whidi can be conveniently and faithfully 
constructed by triangles. In the smooth parts of a surface, 

60 rectangle patches are also commonly used. 

The proposed topological (or connectivity) coding 
scheme described herein was tested on the entire model set. 
The compression performance versus the size of the model 
in terms of the number of vertices is shown in FIG. 12. From 

65 this figure it can be appreciated that, generally speaking, the 
compression performance is better as the size of the model 
is larger. However, since models can have a wide variety of 
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topological structures, the model size is not the only factor 
ctmlrolling the compressioD performance. Statistics about 
the performance of connectivity coding arc summarized in 
Table 1 as follows: 
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As can be seen^ 2.18 bits per vertex on average can be 
achieved. If coding bits are normalized with respect to the 
number of polygon of a mesh there are Z02 bits per polygon, 
For triangular meshes, the average bits per triangle can also 
be computed and is equal to 1.46 bits per triangle. 

The proposed geometry coding scheme herein described 
was tested on the entire test data set. Since the algorithm has 
an embedded property the decoder can rebuild the model at 20 
any bit rates. FIGS. 13 {a) (b) and (c) show the compression 
performance at bit rates of 10, 20, and 30 bits per vertex 
respectively. FIG. 13(i/) illustrates the tradeoff of bit rate to 
distortion. Some statistics based on the entire model set are 
summarized in Table 2: 25 
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FIGS. 14(fl), 14(fc), 14(c) and 14(d) illustrate the resuh of 
constructing the model "dodge" usic^ the different bit rates 
10 bitsA^ertex, 15 bits/vertex, 20 bits/vertex, and 25 bits/ 
vertex, respectively. 

One noticeable artifact associated with high compression 40 
ratios is the "cradc" effect, which refers to gaps between 
different components of a model that are supposed to be 
stitched together. The reason for the crack phenomenon is 
that two vertices £nom different components may have the 
same position in the original model but they are predicted by 45 
different ancestors and, therefore, have different residues 
and get different quantization error. When decompressed, 
this difference caxises the two vertices to take different 
positions. However, the crack does disappear when the 
compression ratio is in the medium and low ranges where a 50 
finer quantization step size is adopted. 

Having thus (Ascribed an embodiment of the invention, it 
be understood by those skilled in the art that many changes 
in construction and widely differing embodiments and appli- 
cations of the invention will suggest themselves without 55 
departing from the spirit and scope of the invention. Thus, 
the disclosure and the description herein are purely illustra- 
tive and are not intended to be in any sense limiting. 

What is claimed is: 

1. A method of encoding a three dimensional model 60 
represented by polyhedral meshes to provide a compressed 
data representation of the three dimensional model, said 
method comprising the steps of: 
A. treatii^ the mesh structure according to the fDllowing 
sequence: 65 
1. select a starting node and a principle link I 
associated with the starting node n^; 
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Z traverse each link associated with the starting node 
in a counter clockwise direction starting with the 
principle link; 

3. for each link traversed in Step 2, add the traversed 
link aiding node n^ to a queue and set the link as the 
principle link, if the eodir^ node is not already in the 
queue; 

4. after each link associated with the starting node n^ 
has been traversed go to the next node in the queue 
and repeat Steps 3 and 4 imtil all nodes in the queue 
have been traversed, and then terminate; 

B. recording each link associated with each node of the 
mesh traversed. 

2. The method of claim 1 further compriang the step of 
identifying the ending node of Step A-3 as a brandi node 
if not in the queue and a merger node if in the queue. 

3. The method of claim 2 further comprising the step of 
identifying each merger node according to a global indexir^ 
scheme. 

4. The method of claim 2 further comprising the step of 
identifying each merger node according to a local indexing 
scheme. 

5. The method of claim 4 wherein the local indexing 
scheme comprises a clockwise tracing or a counter clock- 
wise tracing to the merger node from the starting node n^ 
whichever is shorter. 

6. The method of claim 1 further comprising the step of 
specifying the attributes of the node for each node traversed 
in step A. 

7. The method of claim 6 wherein each of the attributes of 
the node are selected from the group consisting of position, 
normal, color and texture. 

8. The method of claim 6 further comprising the step of 
coding the attributes of the node, comprising the steps of: 

1. ordering the attributes to be coded; 

2. producing a sequence of prediction residues for each 
coded attribute; and 

3. quantizing and coding the prediction residues. 

9. A method of encoding a three dimensional model 
represented by polyhedral meshes to provide a compressed 
data representation of the three dimensional model, said 
method comprising the steps: 

coding the structure data of the polyhedral meshes rep- 
resenting the three dimensional model, wherein the step 
of coding the structure data composes the steps of 
selecting a starting node 11, and a principle link asso- 
ciated with the starting node n_^ and traversing a 
selected link associated with the starting node n ^ in a 
counter clockwise direction starting with the princ^le 
link; 

coding the attribute data of the polyhedral meshes repre- 
senting the three dimensional model; and 

multiplexing the encoded structure data and the encoded 
attribute data into a single data stream. 

10. The method of claim 9 wherein the step of coding the 
stmcture data comprises the step of: 

adding the selected link ending node to a queue and 
setting the link as the principle link if the ending node 
n. is not already in the queue. 

11. The method of claim 10 further comprising the ^ep of 
traversing all nodes in the queue by repeating said traversing 
step and said adding step. 

12. The method of claim 10 further comprising the step of 
idemifying the ending node n^ as a branch node if the ending 
node n^ is not in the queue and as a merger node if the 
identifying node n^ is in the queue. 
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13. The method of claim 12 further comprising the step of 26. The method of claim 25 wherein each of the attributes 
ideatifying eadi merger node according to a gbbal indexing of the node are selected &om the group consisting of 
scheme. position, normal, color and texture. 

14. The method of claim 12 further comprising the step of 27. The method of claim 25 further comprising the step of 
identifying each merger node according to a local indexing 5 coding the attributes of the node, comprising the steps of: 
scheme. ^* ordering the attributes to be coded; 

15. The method of claim 14 wherein the local indexing 2. producing a sequence of prediction residues for each 
sdieme comprises a clockwise tracing to the merger node coded altnbute; and 

from the starting node n,. 3. quantizing and coding the prediction residues. 

16. The method of daim 14 wherein the local indexing lO 28. A method of encodtti^ a three dimeisional model 
scheme comprises a counter clockwise tracing to the merger represented by polyhedral meshes to provide a compressed 
Dode &om the starting node n representation of the three dimensional model, said 

17. The method of daim 9 wherein the step of coding the comprisii^ the steps: 

attribute data comprises the step of specifyii^ the attributes ^ structure data of the polyhedral meshes rep- 

of each node in the polyhedral meshes. 15 resenting the three dimenaonal model, wherem the step 

IS.Tltemethodofclaimnwhereineadioftheattributcs of coding the structure data comprises the steps of 

of each node in the polyhedral meshes are selected from the ^^f^^ iT^t ^"^^^^ ^ 

r vr t 1 J * ^ ciated with the starting node n- and traversing a 

^ta^^°^!?L S;'^'!!^'^ color and texdue^ ^^^^ ^ assodated with the ^g node njk a 

19. T^emethodofckmil7forthercom^ clodcwisc direction starting with the prLiple Knk; 
coding the attributes of each node, compnsmg the steps of: 20 ^^^^^ ^^^^ polyhedral meshes reprc- 

ordcnng the attributes to be coded; senting the three dimensional model; and 

producing a sequence of prediction residues for each multiplexing the encoded stmcturc data and the encoded 

coded attribute; and attribute data into a single data stream, 

quantizing and coding the prediction residues. ^ 29. The method of claim 28 wherein the step of coding the 

20. A method of encoding a tree dimensional model structure data comprises the step of: 

represented by polyhedral meshes to provide a compressed adding the selected link ending node n^ to a queue and 

data representation of the three dimensional model, said setting the link as the prindple link if the ending node 

method comprising the steps of: n, is not already in the queue. 

A. treating the mesh strucmre according to the following 30 30- The method of claim 29 further comprising the step of 
sequence: traversing all nodes in the queue by repeating said traversing 

1. select a starting node n, and a principle link I step and said adding step. 

assodated with the starting node n^; 31. The method of claim 29 further comprising the step of 

2. traverse each link associated with the starting node identifying the ending node n, as a branch node if the ending 
in a clockwise direction starting with the prindple 35 node n^ is not in the queue and as a merger node if the 
tiok- identifying node n^ is in the queue. 

3. for each link traversed in Step 2, add the traversed 32. The method of claim 31 further comprising the step of 
link ending node to a queue and set the link as the ideatifying each merger node according to a global indexing 
prindple hnk, if the ending node is not already in the scheme. 

queue; 40 33, The method of claim 31 further comprising the step of 

4. after each link assodated with the starting node n_, identifying each merger node according to a local indexing 
has been traversed go to the next node in the queue scheme. 

and repeat Steps 3 and 4 until all nodes in the queue 34. The method of claim 33 wherein the local indexing 

have been traversed, and then terminate; scheme comprises a counter dockwise tracing to the merger 

B. recording each link assodated with each node of the 45 °ode from the starting node n^ 

mesh traversed method of claim 33 wherein the local indexing 

21. The method of daim 20 further comprising the step of scheme comprises a clockwise tracing to the merger node 
identifying the ending node n. of Step A-3 as a brandi node ^^f" ^ starting ncxie ^ ^ . , .... 
if not in the queue and a merger node if in the queue. 36. The method of daim 28 wherem the step of codmg the 

22. The method of daim 21 further comprising the step of 50 attribute data comprises the step of specifying the attributes 
identifying each merger node according to a gtobal indexing node m the polyhedral me^es. 

g^jjjgjjjg 37. The method of claim 36 wherein each of the attributes 

23. TTie method of cUim 21 further comprising the step of «*f polyhedral meshes are selected from the 
identifying each merger node according to a local indexing g^up consistmg of position, normal, cotor and texture, 
g^^jjjg 5^ 38. The method of claim 36 fiirther comprising the step of 

24. The method of daim 23 wherein the local indexing coding the attributes of each node, compriar^ the steps of: 
scheme comprises a counter clockwise tracing or a clock- ordering the ^tributes to be coded; 

wise tracing to the merger node from the startic^ node n^ producing a sequence of prediction residues for each 

whichever is shorter. coded attribute; and 

25. The method of claim 20 fiirther comprising the step of 60 quantizing and coding the prediction residues, 
^cifyii^ the attributes of the node for each node traversed 

instepA- ***** 
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